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Relativistic fermions on a lattice are shown to correspond 
to the fluctuations in the localized nonrelativistic fermions. 
Therefore, in contrast to nonrelativistic case, the relativistic 
fermions are critical with universal exponents described by 
the strong coupling limit of the nonrelativistic problem. The 
fluctuations also describe anisotropic spin chain at the onset 
to long range magnetic order whose universality class is the 
Ising model. This generalizes the universality in spin mod- 
els to include multifractal exponents. Finally, analogous to 
the nonrelativistic case, the relativistic fermions may exhibit 
ballistic character due to correlations. 

PACS numbers: 05.45. +b,63. 50. +x,64.60.Ak 



The problem of Dirac fermions with random mass has 
been the subject of various theoretical studies. (^-^] For 
weak disorder, the model describes one-dimensional (ID) 
spinless fermions with random hopping matrix elements 
]2|, a toy model that exhibits localization transition in 
ID. Both the spin-Peierls and the spin-ladder systems 
with non-magnetic impurities have been described by 
Dirac equation with random mass. (^^] The zero-energy 
states which dominate the low temperature thermody- 
namic have been the subject of main focus as it provides 
an exactly solvable model whose multifractal spectrum 
can be computed analytically . ^,|| It turns out that in 
contrast to the usual ID localization problem for non- 
relativistic fermions, correlation functions for relativistic 
fermions exhibit power-law decay with universal expo- 
nents, n 

The central result of this paper is that the differences 
between the relativistic and the nonrelativistic fermions 
in presence of disorder, can be attributed to simple rela- 
tionships between these two problems. We show that rel- 
ativistic fermions are described by the fluctuations in the 
exponential decaying wave functions of the nonrelativistic 
fermions. New insight into both the absence of localiza- 
tion and the universality of zero mode Dirac fermion is 
gained by their mapping to the anisotropic spin- 1/2 chain 
at the onset to long range order(LRO). Zero mode Dirac 
fermions describe critical spin chains with O(l) symme- 
try while the corresponding non-relativistic fermions are 
related to 0(2) spin chains. We argue that the massless 
excitations of the spin chain responsible for long range 
magnetic correlations provide mechanism for derealiza- 
tion of relativistic fermions and the statement of univer- 
sality is the well known statement that the anisotropic 
spin chains belong to the universality class of the Ising 
model. Although the setting we describe is quite general 



in the context of disordered systems, for concreteness we 
will consider quasiperiodic disorder which exhibits local- 
ization transition in nonrelativistic lattice problem. The 
reason for this particular choice is that the results from 
exact renormalization group (RG) § as well as rigor- 
ous mathematical analysis H can be used to compare 
and contrast the role of disorder in relativistic vs non- 
relativistic problem. For quasiperiodic disorder, non- 
relativistic problem known as Harper equation fiof has 
been studied extensively. Another interesting result of 
this paper is that it provides a new way to understand 
recently discussed strong coupling fixed point of Harper 
equation which describes the universal aspects of the lo- 
calized phase. The strong coupling universality is the 
statement that the universality class of anisotropic spin 
chain is determined by the Ising limit. 

We consider a nearest-neighbor (nn) tight binding 
model (tbm) of spinless non- interacting fermions in a dis- 
ordered potential V„, 



H = ^[<4 +1 c„ + 2XV n c! n c rt 



(1) 



Here the c's are canonical Fermion operators and A is the 
strength of the aperiodic potential V(x). The eigenstates 
of this lattice system are given by the following discrete 
Schrodinger equation, 



lpn-1 + i'n+l + 2AVvVn = Elj] ri 



(2) 



For disordered systems exhibiting exponential localiza- 
tion, the fermion wave functions ip n can be written as 



1pn 



(3) 



where 7 is the inverse localization length 7 = £ . The 
tbm describing the fluctuations rj n in the exponentially 
decaying envelope is given by the following pair of equa- 
tions, 



e 'Vn+1 + e 'V„-l 

e 7 ^ +1 + e" 7 ^.! 



2AK< t 



Er) r n 
Er, l n 



(4) 



Here rf and rf respectively describe the fluctuations 
to the right and to the left of the localization center, 
which is chosen to be at n = 0. As shown below, these 
equations describe relativistic fermions on a lattice for 
zero energy states as in the long wave length limit, the 
equations reduce to the Dirac equation. We replace n 
by x and write rj n ±x = e :iz ' lp i](x), where p is the momen- 
tum canonically conjugate to x. The equation (4) for the 
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fluctuations for E = state can be described by the fol- 
lowing non-Hermitian Hamiltonian Hfi uc and its adjoint, 

H f i uc = e~^e lp + e 7 e~ 4p + 2XV(x) (5) 

In the limit (p — ► 0) , the system for £ = reduces to the 
Dirac equation, 

[ga x p - i(2XV(x) + 2cosh(~f))a y }ri{x) = (6) 

where f](x) is a two-dimensional spinor r)(x) — 
(rf(x), rf{x)). Therefore, the two-component structure 
of Dirac spinor arises naturally when we consider fluc- 
tuations about exponentially localized wave functions. 
Here g = 2sinh(p/) is the velocity of the Dirac fermions 
while the mass of the Dirac fermions m(x) — 2((\V(x) + 
cosh("f)). The o^, k = x,y are the Pauli matrices. 

Therefore, on a lattice, the Dirac fermions with dis- 
ordered mass are the fluctuations of the nonrelativistic 
localized fermions. This would imply the absence of ex- 
ponential localization for relativistic fermions provided 
the equation (4) has a solution with E = 0. 

Before we discuss the solution of the discrete Dirac 
equation, we show the equivalence between the disor- 
dered Dirac Hamiltonian and a spin chain with disor- 
dered magnetic field. The anisotropic XY spin- 1/2 chain 
in a transverse magnetic field is, 

H = -H^ a «+^ + ^M+i + 2AK<]. (7) 

The e~ a and e a respectively describe the exchange in- 
teractions along the x and the y directions in spin space 
and therefore, the spin space anisotropy g is given by 
g = 2sinh(a). It is well known that the Jordan- Wigner 
transformation, |]l6f transforms the spin problem to spin- 
less fermion problem where fermions are the quasiparticlc 
excitations of the spin chain. The eigenstates of the exci- 
tations are described by the following coupled equations, 

e~ a Vn+i + e a ti-i + 2XV„v} = Eif n 

e a vl+i + e-Wn-x + 2\V n r,l = Eri (8) 

For E = 0, these equations reduce to equation (4) with 
the parameter a equals the inverse localization length 7. 
However, E = is the solution of this equation only 
at the onset to LRO, and therefore, for anisotropic spin 
chain with disordered field V n , this mapping between the 
spin and fermions provides an interesting way to describe 
the onset to LRO: the massless excitations of the critical 
anisotropic spin chain describes the fluctuations in the 
exponentially localized excitations of the isotropic chain 
provided the anisotropy parameter is the inverse local- 
ization length. Alternatively, the mapping provides a 
new method to determine the localization length of tbms 
in the presence of disorder. For the case of quasiperiodic 



disorder, V n = cos(27r(ern + 0)), ( where a is an irrational 
number and is a constant phase factor), the localiza- 
tion length is = log(A). fTTi] Therefore a = log(X) 
is the equation for the critical line, a result known from 
earlier numerical study |ll]] for the spin chain. 

The Ising limit (a — > 00) of the critical spin chain 
coincides with the A — » 00 limit of the Harper equa- 
tion. This is the strong coupling limit of the Harper 
equation ]tJ that has been studied by exact RG as well 
as more rigorous analysis The existence of a strong 
coupling fixed point shows the universality of the fluc- 
tuations in the exponentially decaying solutions in the 
localized phase. The equivalence between the strong cou- 
pling limit of Harper equation and the critical Ising spin 
in a quasiperiodic transverse field provides a new way to 
understand this fixed point: the statement of universal- 
ity of the strong coupling fixed point of Harper equation 
is equivalent to the statement that the anisotropic spin 
chain in a quasiperiodic field belongs to the universality 
class of the Ising model, a result well known for periodic 
spin systems in arbitrary dimension. 

We now focus on the E = solution of the Harper 
equation in the localized phase as the fluctuations about 
this describe the Dirac fermions. The universal features 
of this solution can be studied by an exact RG scheme. @] 
The earlier studies of this problem was confined to the 
band edges of the eigenspectum. To show the univer- 
sal features of the solution, with a equals inverse golden 
mean, all sites except those labeled by the Fibonacci 
numbers F m are decimated. At the m th decimation level, 
the tbm describing the fluctuations is expressed in the 
form 

(n) Vn . (9) 

The additive property of the Fibonacci numbers provides 
an exact recursion relations for the decimation functions 
e m and f m . In the strong coupling limit f m is found to 
approach zero simplifying the RG flow as, W 

e m +i(n) = -e m _i(n + F m )e m (n) (10) 

Numerical iteration of this equation provides an ex- 
tremely accurate method to distinguish extended, lo- 
calized and critical states. The extended and localized 
phase correspond to trivial asymptotic RG flow while 
the critical behavior is found to be characterized by a 
nontrivial asymptotic 6-cycle in the decimation function: 

Figure 1 shows the self-similar fluctuations about ex- 
ponentially localized solutions in Harper equation, which 
describes zero-mode Dirac problem with quasiperiodic 
mass. In the strong coupling limit, which also describes 
the massless mode of the Ising chain the wave function 
is given by the following algebraic equation, 

n 

5n ^ exp-Mn)} = exp-[T,log(2\Vi\)] (11) 
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For random disorder, the above summation can be done 
explicitly resulting in exact solution [||J|] , that facilitated 
an exact calculation of multi-fractal spectrum. 

In order to seek the physical meaning of universal 3- 
cycle, and see whether the statement that the critical 
exponents of spin chains are determined by the symmetry 
of the spin chain applies to the multifractal exponents, 
we compute the f(a) curve describing the multifractal 
spectrum associated with the self-similar wave function 
or the inverse participation ratios P, 



P(q,N) 



E\Vn\ 2 

_ dr 
dq 

f(a) = aq- r(q) 



N— T <.<i) 



(12) 



The free energy function r(q) and its Legendre transform 
f(a) were found to be A independent only for for positive 
values of q and hence only left half of the f(a) curve is 
universal. Therefore, for quasiperiodic spin chains at the 
onset to LRO, scaling exponents for only positive mo- 
ments of the participation ratio are determined by the 
spin space symmetry. 

We next show that the presence of correlated disor- 
der can result in propagating solutions for the disor- 
dered Dirac equation, as is known to be the case for the 
non-relativistic fermions. [ p"3"| For quasiperiodic disorder, 
dimer-type correlations can be introduced by replacing 
9 n = 2nan in V(0 n ) by the iterates of the supercritical 
standard map, describing Hamiltonian systems with two 
degrees of freedom. |U5lE3 



Jn+l 



28 L 



K 

-— sin(27r6>„). 
Zn 



(13) 



We use iterates that describe golden-mean cantous, 
which has been shown to exhibit dimer-type correlations, 
and leads to Bloch-type states for the nonrelativistic 
fermions. [|l4| Here, we will confine to the Ising model, 
described by, 



Vn+l 



2\cos{2n6 n )r 1 l n = Erfc 
_! + 2Xcos(2Tr9 n )r 1 l = Erfc 



(14) 



We determine the critical A as a function of K, and ana- 
lyze the massless mode of the Ising model using the RG 
methodology 0. As shown in figure 3, nontrivial 6-cycle 
degenerates to trivial fixed points at certain special pa- 
rameter values. The origin of these trivial fixed points, 
has been traced to a hidden dimer in the quasiperi- 
odic iterates describing the golden-cantorus. H] At these 
points, the relativistic massless mode of the Ising model 
is ballistic. Further details of the effects of correlated 
disorder on magnetic and spectral transitions will be dis- 
cussed elsewhere. 



We would like to point out that although the rela- 
tionship between the fermion and the spin problem is 
well known, what is new here is the that disordered spin 
chains provide important key in understanding the role 
of disorder in transport properties of relativistic and non- 
relativistic fermions. The differences between these two 
fermions is traced to the differences between spin chains 
with 0(1) and 0(2) symmetry where the spin chains 
with O(l) symmetry exhibits magnetic long range cor- 
relations. The correspondence with the spin chain also 
explains why E = mode of Dirac equation is special: 
the root of criticality of E = Dirac fermions can be 
attributed to the fact that it corresponds to infinite mag- 
netic correlation length. We argue that these magnetic 
correlations account for the absence of exponential local- 
izations in relativistic fermions. It is rather interesting 
that the relativistic fermions can be viewed as the fluc- 
tuations about localized nonrelativistic fermions and this 
picture also explains the absence of localization of Dirac 
fermions. Furthermore, the universal aspects of relativis- 
tic fermions are described by the strong coupling limit of 
the nonrelativistic fermions. This universality implying 
that the positive moments of the participation ratio have 
universal exponents , generalizes the well known state- 
ments of universality of periodic spin chains to quasiperi- 
odic spins exhibiting multifractal exponents. Finally, we 
would like to mention that many of our conclusions are 
valid for random as well as pseudorandom disorder. p7| 
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FIG. 1. (a) Absolute value of the fluctuations for Harper 
equation for E — states with (j> = .25 . At the Fibonacci 
sites, we see period-6 behavior(period-3 in absolute value): 
VlFm \ = l^fWal wmcn ls > independent of A. 



FIG. 2. Numerically obtained (a) f(ct) curves for A — > oo( 
solid curve) and A = 1.5 (lines with crosses). 

FIG. 3. (a)Critical A as a function of K for the Ising model, 
(b) RG 6-cycle showing the variation in the renormalized cou- 
pling e after Fibonacci decimation. 
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